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Abstract 

We investigate the collision-limited electronic Raman response and the attenuation of 
ultrasound in spin-singlet d-wave superconductors at low temperatures. The dominating 
elastic collisions are treated within a t-matrix approximation, which combines the de- 
scription of weak (Born) and strong (unitary) impurity scattering. In the long wavelength 
limit a two-fluid description of both response and transport emerges. Collisions are here 
seen to exchisivcly dominate the relaxational dynamics of the (Bogohubov) quasiparticle 
system and the analysis allows for a clear connection of response and transport phenomena. 
When applied to quasi-2-d superconductors like the cuprates, it turns out that the 
transport parameter associated with the Raman scattering intensity for Big and B2g photon 
polarization is closely related to the corresponding components of the shear viscosity tensor, 
which dominates the attenuation of ultrasound. At low temperatures we present analytic 
solutions of the transport equations, resulting in a non-power-law behavior of the transport 
parameters on temperature. 

PACS: 67.57.Hi 74.20.-z 74.20.Fg 74.20.Rp 74.25.Fy 74.25.Ld 

1 Introduction 

During the last few decades a large variety of so-called unconventional superconductors 
have been discovered, among these the superfluid phases of ^He [1-4], the heavy Fermion 
systems [5-7], the cuprates [8, 9] and the Ruddlesden-Popper system Sr2Ru04 [10]. The 
unconventional pairing correlations in these systems manifest themselves on the one hand in 
an anisotropy of the pair potential or energy gap of less symmetry than the underlying band 
structure or in the occurrence of additional spontaneously broken symmetries besides the U{1) 
gauge symmetry. On the other hand their existence can be detected from the sensitivity of 
thermodynamic quantities like the transition temperature or the equilibrium energy gap to 



even small amounts of non-magnetic impurities. 

In an earlier publication [11], a simple two-fluid description for these unconventional su- 
perconductors was formulated, which emerges from the BCS theory in the long wavelength 
(q — > 0) and stationary (a; — > 0) limit, sometimes also referred to as the local equilibrium. 
As a result, even analytical results were obtained for the temperature dependence of the 
local reactive response functions of the normal component, the Bogoliubov quasiparticles 
(specific heat capacity, spin susceptibility) and the condensate (superfluid density, magnetic 
penetration depth) . Clearly, a comprehensive two-fluid description of superconductors should 
describe the more general situation beyond local equilibrium and should therefore contain 
the dissipative response of both the quasiparticle system and the condensate. A flrst step in 
this direction was the derivation of a general form of a certain class of quasiparticle transport 
parameters for unconventional superconductors in reference [12]. The results for the impurity- 
limited transport of momentum (shear viscosity) and energy (diffusive thermal conductivity) 
were, however, discussed exclusively for the case of superfluid ^He-B with silica aerogel 
forming the impurity system. Moreover, the two aspects of response and transport tradition- 
ally appear as fairly remote aspects of the reactive and dissipative dynamics of a superconductor. 

Therefore, this paper is devoted to a unified description of response and transport in su- 
perconductors. To be specific, we limit our considerations to unconventional spin-singlet 
superconductors with d-wave pairing correlations, in view of an applicability to hole-doped 
cuprate superconductors. We would furthermore like to concentrate on the electronic Raman [13] 
and stress tensor response. When treated in the long wavelength limit, the corresponding 
response functions can be shown to be separable into normal (Bogoliubov quasiparticles) and 
superconducting (pair condensate) contributions, respectively, hence allowing for a two-fluid 
description at arbitrary quasiclassical frequencies in the homogeneous (q — > 0) limit. While 
the dynamics of the condensate can be characterized by some pseudo-conservation laws 
governing the macroscopic phase of the order parameter (reactive response), as well as by 
pair breaking processes (dissipative response), the system of Bogoliubov quasiparticles shows 
purely relaxational behavior in the long wavelength limit. For a quantitative study of the latter 
behavior we consider impurity-limited transport, believed to dominate at low temperatures. 
Collisions of the BQP with impurities are treated within the t-matrix approximation [14,17-19], 
which is limited here to s-wave scattering. The description thus allows one to treat the cases 
of weak scattering (Born limit, 5o 0) and strong scattering (unitary limit Sq 7r/2) on 
the same footing. The theory is applied to an analysis of the electronic Raman response and 
the attenuation of ultrasound in d-wave superconductors. Our formulation is general enough 
to include the aspects of universal transport, which was discussed previously in the literature 
in context with electronic conductivity [21], diffusive thermal conductivity [22], ultrasound 
attenuation [23] and electronic Raman response [24] . 

The paper is organized as follows: After a discussion of the equilibrium properties of unconven- 
tional superconductors in section 2, we establish a general response theory in section 3 which 
is based on the classiflcation of external perturbation potentials through vertex functions Cp, 
which correspond to the specific experiment, testing the response. Section 4 then deals with 
the derivation of a two-fiuid description of the response in the long wavelength limit (q 0) 
at arbitrary quasiclassical frequencies u. The effects of the long-range Coulomb interaction are 
explicitly taken into account in the derivation of the response functions. Results for the con- 
densate response, which are known from the literature, are briefly rederived for completeness 
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at the end of this section. Section 5 is devoted to the response and relaxation properties of the 
system of thermal excitations of the superconductor, the Bogoliubov quasiparticles (BQP). It is 
shown that the dynamics of the macroscopic density fluctuations of the BQP system is entirely 
relaxational and describable by a set of macroscopic relaxation times, which depend on the ver- 
tex function Op. This concept allows for the derivation of quite general equations, which relate 
the response functions of the BQP system to the corresponding transport parameter of given 
vertex function dp. The influence of the long range Coulomb interaction on the qualitative form 
of the transport parameter is studied. In section 6 we consider the special case Op = 1, i. e. the 
relaxtion of the macroscopic quasiparticle density in the absence of the Coulomb renormalization 
and make contact to earlier work on intrinsic density relaxation and the second viscosity. In 
Section 7 we derive the explicit form of the impurity-limited relaxation time, which enters the 
transport parameters using the t-matrix approximation for the impurity self-energy. In sections 
8 and 9 we discuss the similarities in the temperature dependence of the transport parameters 
associated with electronic Raman scattering and ultrasound attenuation for various parameters 
characterizing impurities in the weak (Born) and strong (unitary) scattering limit. Section 10 
is finally devoted to our summary and conclusion. 



2 Equilibrium properties 

It is well established that the pairing correlations in cuprate superconductors are unconventional 
in the sense that the Fermi surface average of the gap function Ap vanishes, i. e. 

(Ap)FS = (1) 

As a special form of the gap anisotropy we consider the case of Big gap symmetry [9], 

Ap = Ao(T) cos(2</.) (2) 

Note that the nodal structure of such a gap function implies, that the thermal excitations of the 
system, the Bogoliubov quasiparticles (BQP), which have an excitation spectrum of the usual 
form 

Sp = y^i + A2 , (3) 

can be created at arbitrary small energies (nodal quasiparticles). This is reflected in the 
form of the BQP density of states Ns{Ep) 

Ns{Ep) _2^^fAo\ i?p-.o Ep 



No TT \EpJ Ao ^ ^ 

varying linearly in the quasiparticle energy at low energy in the clean limit. In Eq. (4) K 
denotes the complete elliptic integral of first kind and Nq is the electronic density of states at 
the Fermi surface for one spin projection. The statistical properties of the excitation gas can 
conveniently be described by the thermal Fermi-Dirac distribution 

(5) 



Note that a chemical potential term is missing from (5) , since the number of thermal excitations 
is not fixed. 
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3 Response theory 



In order to test the response and transport properties of a superconducting system, one has to 
apply external perturbation potentials, which can be classified in the following way: 

^'^ij = X] "k<^^a (6) 
a 

The first and second term in (6) describe the coupling of the electronic system to the 
electromagnetic scalar ($) and vector (A) potential, respectively. The third term represents 
the electronic coupling to a typical Raman scattering process, with an incoming photon of 
energy JIloj, momentum ^k/ and polarization along A^, and a scattered photon with energy 
hujs, momentum ^k^ and polarization along A*^ leaving the sample. This process couples, for 
example, to electronic excitations near the Fermi surface with energy transfer hu = fkoj — Hus 
and momentum transfer = ^k/ — ^kg, and is describable by a k-dependent so-called 
Raman tensor 7i^., which we have approximated in Eq. (6) by the inverse effective mass tensor 
{Mj^^j^jy = d'^^]f_/lT?dk^dk,y. In (6) rg = /mc? denotes the Thompson radius. The fourth 
term in (6) represents the electronic coupling to the lattice strain field 5uij, which leads to a 
dissipative response of the electronic stress tensor and hence to the attenuation of ultrasound. 
In (6) Vk = dE]f^/hd\s. = (^k;/-E'k)vk is the group velocity of the BQP and d the dimension of the 
system. The r.h.s. of Eq. (6) generalizes the k-space structure of the perturbation potentials 
by introducing a k-dependent so-called vertex function Ok, together with a collection of fictive 
potentials 5^a, related to each vertex. In this spirit, the vertex function Ok is related to the spe- 
cific experiment (electromagnetic response, Raman response and relaxation, sound attenuation) 
under consideration. In the Raman case one has flk = 7k = e/ -Tk "^S with ej^s the unit vectors 
in the direction of A^'*^ and 6^^ = rolA^HA"^!. In the case of sound attenuation one may write 
cik = q • fk • u with the definitions piV]^j - p • \'k6ij/d = PFVpi^k/ Ekjcr^ij = {^.k/Ekja^ij and 
^^(7 = |q||u| [25,26]. It is interesting to note, that in quasi-2-d systems the vertex functions 
Ok for Raman scattering and sound attenuation coincide in case of Bi^- (ak = cos(2(i))) 
and B2g- (ak = s'm(2(p)) symmetry. As we shall demonstrate explicitly in section 9, these 
Raman polarizations can be shown to correspond to the attenuation of transverse sound, i. e. 
q _L u, if q is oriented parallel (B2g-symmetry) to the crystal (a-) axis, or is tilted by 45° 
(Bi^-symmetry) from it. 

In the homogeneous limit q ^ the response and transport properties can be clearly separated 
into a condensate and a BQP contribution. In the case of electronic Raman scattering, the 
condensate contributes to what is referred to as the pair-breaking Raman effect, which has been 
extensively discussed in the literature [27,28]. In contrast, the condensate does not contribute 
to dissipative processes like, for example, the impurity-limited Raman effect or the attenuation 
of ultrasound. 

The total density fluctuation of the superconductor can, as usual, be written as: 

5nk(q, u;) = 5 (4+qaCka) (uj) (7) 
The presence of the perturbation potentials (6) gives rise to a macroscopic density response 

Sna{q.,oj) = ^^ap5np{ci,Lv) = {ap5np{q_,u})) (8) 

pa 

(...) = ^E... w 

pa 



e$ - -Vk • A + m{M^^}ijroAlA 



PiVkj 



pVk 
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In the presence of the long range Coulomb interaction 

K, = ^ (10) 

the density fluctuations (5n(q, w) = (5ni(q, a;) give rise to a molecular potential 

,56 = Vq<5ni(q,a;) (11) 
which adds to the external potentials 6^^^ and leads to the coupled response 

5na(q,a;) = xi°i(q, ^)<5ea + xi?(q, c^)<56 (12) 
Sni{ci,uj) = x£'(q,a;)(5Ca + Xn (q,'x;)(56 

The system of Eqs. (12) can easily be solved with the result 

5na{q.,aj) = Xaa(q, t^)'5Ca 

xl?^q,w) 



Xaa(q,u;) = Xaa(q,w) (q) . 

Xii (q,^^) 

e(q,w) = l-yqXn(q,a;) 



e(q,a;)_ 



(13) 



4 Two— fluid description 

It can be shown that a two fluid description emerges close to the long wavelength limit q ^ 0, 
i. e. (5nk(q —>■ 0,a;) can be decomposed into a condensate (Sn^) and a quasiparticle (Sn^) 
contribution (assumig at this stage, that the collisions are not pair-breaking): 

5nk(q -> 0, w) = 6nl{uj) + Sn^{u;) (14) 

The condensate contribution (pair response) has the gauge-invariant form [29] : 

6nliuj) = -Ak(a;) l^d^l^^ - '-hw6^^ (15) 

with Ak(a;) the Tsuneto function in the long wavelength limit 

4A20k . 1 Ek 
= 4lF-^ ' = 2^ W ^''^ 

and ,5^^^^ denotes the part of (6) with Ok = a_k- In Eq- (15) Sip denotes the nonequilibrium 
phase change of the order parameter, the time derivative of which is connected with the external 
perturbation potentials of even parity through the Hamilton-Jacobi (or generalized Josephson) 
relation 

-.6, ^-,^^6^= = + E (17) 

where the short-hand notation (9) for the momentum sums has been used. The physical inter- 
pretation of Eq. (17) as a Josephson relation has been emphasized by writing out explicitly the 
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contribution from the scalar potential $ to the phase change. The other terms in (17) will turn 
out to vanish except for the Ai^ Raman polarization, to be discussed later. Inserting (17) into 
(15) leaves us with 



^'nl_ici^O,uj) = -Ak < 



Note that the pair response vanishes for a constant vertex Op = const as a consequence of 
the gauge invariance of the theory, expressed through the relation (17). The total generalized 
response function Sua can be decomposed into its pair (P) and quasiparticle (Q) contributions 
as follows 

Sua = Sn^ + Sn^ 
5n^'Q = {apSnl'^) = x^ak'^S^a (19) 

Xaa Xaa ~^ Xaa 

^From (18) the renormalized pair response function Xaa be written in the form 

xL = -Aaa+^ (20) 
All 

^ab = (ApCp^p) 

It was shown in rcf. [29] that in the absence of collisions the quantity Qx^ai'^) entirely describes 
the Raman response of the superconductor, the so-called pair-breaking Raman effect. 



5 Homogeneous quasiparticle transport and relaxation 

In what follows we shall therefore concentrate on the BQP contribution (normal component 
in the spirit of a two-fluid description) to the response, transport and relaxation properties. 
Restricting our consideration to the case of Raman scattering and sound attenuation, the vertex 
Ok has positive parity, i. e. a_k = Ok- In this case, using Eq. (9), one may define a macroscopic 
BQP density via 

5n^ = (^ap^Si.p'j , (21) 

with (5fp the deviation of the BQP distribution function from equilibrium. Special cases include 
then the BQP Raman response function Srij and the stress tensor response function q • 11*^ • u. 
Here 11*^ denotes the BQP stress tensor and the unit vectors in the direction of propagation (q) 

and polarization (u) emerge from the standard representation of the strain tensor Suij oc QiUj 
[26]. It should be emphasized that the quasiparticle stress tensor is defined as the momentum 

current 

hp = Svp + ypSEp ; yy^ = --^'^^ 



SEp = ^S^f 
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and, strictly speaking, differs therefore from a generalized density. However, we are able to 
show in the appendix, that whereas the reactive response of densities and currents is indeed 
qualitatively different, their dissipative response, i. e. their transport parameters are the same. 
In the long wavelength limit q — > 0, dvp obeys the scalar kinetic equation [30] 



(22) 



where SI^ represents the collision integral for the quasiparticle system. Following ref. [32], we 
decompose the collision integral (5/k into contributions originating from elastic (e) and inelastic 
(i) scattering processes: 



dit 



511 



'p 



2 
P 

p ^p 



(23) 



In Eq. (23) we have used approximate forms for the collision integrals, applicable for distri- 
bution functions of positive parity Sv^i^ = 5fk, which guarantee the Bogoliubov quasiparticle 
number conservation for elastic scattering and allows for describing the fact that the number of 
Bogoliubov quasiparticles is not conserved in context with inelastic scattering processes. In (23) 
Tp ' denote the impurity-limited and the inelastic quasiparticle relaxation times, respectively, 
of the superconductor, the first of which will be specified in more detail in section 7. Before 
we perform the Coulomb renormalization, dictated by Eq. (13), it is instructive to study the 
relevant response functions Xaa^^ for purely elastic scattering 



Xai°^('^) elastic 



(24) 



e 



and purely inelastic scattering 




(25) 



Xr^(c^)inela.tic = -S' „(a;) + ^ 

KbH = 

tq{u)) = 



It is important to note that Eq. (24) expresses the number conservation law Xn'^''(^')clastic = 
for elastic scattering processes in the long wavelength limit. For inelastic scattering, however, 
as represented by Eq. (25), there occurs the well-known phenomenon of intrinsic quasiparticle 
relaxation [30,31], described by the lifetime tq{lo), which is finite below as a consequence of 
the nonconservation of the BQP number density, explicitly built into the inelastic part of the 
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collision integral (23). It is seen to diverge in the limit Ap — since the number of quasipar- 
ticles is conserved in these processes in the normal state. Therefore Eq. (25) is reminiscent of 
a viscoelastic description of the generalized response, which interpolates between the hydrody- 
namic {cotq — > 0) and the coUisionless {lotq — > oo) limit. If both elastic and inelastic scattering 
processes occur simultaneously, the situation becomes more complicated, since one has to solve 
the integral equation 



(+) 



1 — iooT^ ' 1 — ioJT^ 



^P 'P 



+ 



(26) 



Here the index (+) denotes the positive parity of the distribution functions (5z^k a-^d /i^ with 
respect to the operation k — > — k. ^From Eq. (26) one immediately observes, that the relaxation 
rates 



■pe,i 
^ k 



do not simply add up 



= r^k + TL 



(27) 



to result in an effective relaxation time 



^k 



^ k 



r^k + 



(28) 



but there appear mixing terms originating from the collision operator in (26). It should be 
noted that Eq. (26) is a straightforward generalization of the result (3) of ref. [32] to the 
superconducting case. Since in what follows, we are only interested in the homogeneous limit 
q — > of the quasiparticle response, we may follow the argumentation of ref. [32] in solving Eq. 
(26) to get the final result for the full response function (c. f. Eq. (13)) after the Coulomb 
renormalization: 




p - ^^Tp 



(29) 



Note that the terms oc e which describe the complicated mixing of clastic and inelastic 



oo. 



contributions (c. f. ref. [32]), can be neglected in the long wavelength limit q — > 0,e 
Therefore one may state that in the Hmit q — ^ the result for the response function Xaa, 
which includes the effects of the long range Coulomb interaction, has a form characteristic of a 
quasiparticle number conservation law for the BQP, with entering as the effective relaxation 
time, and the phenomenon of intrinsic quasiparticle relaxation becomes more or less irrelevant 
for charged systems except for special experimantal situations described in chapter 5.3 of ref. [31]. 



Having established the response functions of the quasiparticle system, we would next like to 
clarify an important physical consequence of the relaxation equation (22). It turns out that 
from (22) one may derive a set of homogeneous relaxation equations for the BQP densities Sria 
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by multiplying (22) with ap(^p/E'p) and summing on momentum p and spin a. As a result we 
find that the dynamics of the BQP system is entirely governed by relaxation processes on time 
scales set forth by vertex-dependent BQP relaxation times Taa- Assuming the fictive potentials 
5^a to vary as oc exp(iq • r — iuit), these relaxation processes obey a set of general equations for 
each vertex Ok [33] 



r£(a;) 



loc 



which describes the relaxation of the BQP density back to its local equilibrium value 

loc 



- xS(- = 0) = -{(a^-a^) 



(30) 



(31) 



El 



(32) 



Prom Eq. (30) we immediately get 



Xaai^) = 



Q 

Xaa 



(0) 



1 - iuJTaa{uj) 



(33) 



The effective quasiparticle relaxation times r^(a;) are obtained as [33] 



1 



(34) 



With these results, the response function x2i(^^) can be decomposed into its real and imaginary 
parts as follows: 



xZ{^) 



'(a;) 



(35) 



Xaa 



Q 
Xaa 



(o)r; 



Q/ 

aa 



[l + a;r£"]2 + [a;r£']2 

where the prime (/) and the double-prime (//) refer to the real and imaginary part, respectively. 
Eqs. (33, 35) describe the general connection between response and transport of the BQP system 
in the homogeneous limit for a given vertex at- Note that the quantity Taa can be interpreted 
as the generalized quasiparticle transport parameter of the superconductor, since one may write 

^n'i = x^(c^ + I^H^ (36) 

response transport 

in which 5^a and fa = —icoS^a play the role of fictive potentials and forces, respectively. In 
the hydrodynamic limit (w — ^ 0) we obtain the following result for the effective quasiparticle 
relaxation times Taa [33] : 



— 1 2 ^ p * 

a] -^VpTp 



'^PE^Vp 

^Up 



(37) 
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and the transport parameter Taa [33] : 

li- ^a^a(^) = ^(ap - -af ^ypr;^ = ^(ap - af (38) 

Note that the effects of the relaxation time tq (c. f. Eq. (25)), originating from the quasiparticle 
nonconservation in the inelastic scattering channel, have completely disappeared from the result 
for Taa{u} = 0) in the long wavelength limit as a consequence of the long-range Coulomb inter- 
action. The physical consequences of intrinsic quasiparticle relaxation can be best studied for 
neutral pair-correlated Fermi systems, to which we would like to devote the following section. 



6 Intrinsic quasiparticle relcixation and second viscosity 

Although physically relevant only for neutral and not for charged systems, it is interesting to 
investigate the response functions Xaa^^ and the transport parameters Taa^^ in the absence of 
the long-range Coulomb interaction, i. e. before the renormalization manifested through Eqs. 
(12). Restricting the considerations of this section to the clean case (i. e. purely inelastic 
scattering) , the response function is then given by Eq. (25) and reads in the limit a; — (i. e. 
in local equilibrium): 



X, 



Q(0)loc Q(0) 



xr(-=o) 




(39) 



The BQP relaxation equation is of the form (30) 



Sn9 = - 



' aa 



Sn9 - 5n9 



but with a different effective relaxation time 

1 



Xaa 



(0) 



Xaa [i^) 



(40) 



(41) 



entering the representation 



vQ(o) 

Aaa 



■ Q(o)/ X 



(42) 



Eq. (40) represents a straightforward generalization of the problem of intrinsic BQP density 
relaxation occurring in neutral Fermi superfluids [30], which can be obtained from (40) in the 
special case Op = 1. The decay of the BQP density dnf occurs then as a consequence of the fact 
that the number of Bogoliubov quasiparticles is not a conserved quantity in inelastic scattering 
processes. Since the quasiparticle number is conserved in the normal state, the BQP density 
relaxation time must diverge in the limit as T — > T~ [30] . In the hydrodynamic limit {u — 0) 
we obtain results for the effective quasiparticle relaxation times r^^°^ 



limrW)(a;) 



2 

(^pE^yp 



P 



2 ^p 

(^PE^yp 



P Vp 



(43) 
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and the transport parameter Taa^^ '■ 




Urn r£(0) (c) = ( al^y^rl, ) + / (44) 



which are quantitatively entirely different from the corresponding Eqs. (37) and (38). The 
transport coefficient T^^^^^ = (^3 turns out to be nothing but the second viscosity (s, which was 
evaluated for neutral Fermi superfiuids by Wolfie and Einzel [30] : 



It diverges in the limit Ap — >^ as a consequence of quasiparticle number conservation in the 
normal state. When applied to superfiuid ^He-B, the inelastic relaxation time can be taken 

from the appendix 1 of ref. [15]. It should furthermore be noted, that Eq. (26) may serve and has 
already been used as a starting point for a calculation of several relevant transport parameters 
of dirty Fermi superfiuids like ^He in aerogel [12, 16]. 



7 The resonant impurity scattering model 

The aim of this section is to determine the relaxation time Tp of the superconductor in the 
limit of low temperatures where the impurities play the major role for quasiparticle scattering 
processes, i. e. rj^ — > r^. In that case, the relevant scattering parameters are the impurity 
concentration rii and the scattering phase shift Sq (which we would like to restrict to the case 
of s-wave scattering for simplicity) , giving rise to a normal state scattering rate 

^ 2ni . 2. 

— = — — sm So (46) 

In the presence of impurities, the BQP energy gets renormalized through the impurity self- 
energy Se via 

Ep = Ep + Ee(4) (47) 
with Se evaluated within the t-matrix approximation [17] in its self-consistent version [18]: 

Ee(4) = ^ ^(^L_ (48) 

Here, the complex function 

extends the density of states of the superconductor N^{Ey,) /Nq = lkD{Ep) to include impurity 
effects. From the impurity self energy we obtain the elastic scattering rate [34] 



- = ^SS (£ ) = ^^(-^P) cos^(5o + sin2^o|D(^p)p 

^ TN |cOs2 5o + sin2(5oL'2(^p)|2 
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The explicit form for l/Tp becomes particularly simple at high energies and in the limit Ep — ^ 0, 
where the full self-consistent treatment of the renormalization £'k is necessary (Sq = S>Se(0)) 
[20,21]: 



E P 



V"2 



[Si'2+A2] 



3/2 



Ep > S'o' 

Er.^0 



(51) 



8 Impurity— limited transport 



In this section we shall exploit the general result for the transport parameter Taa, derived in 
section 5 further and investigate its dependence on temperature and the parameters describing 
the impurity scattering. Let us recall that 



]^mirf» = r£(r) = l^(ap-a) 

pa 



4 
El 



(52) 



In what follows, wc wish to restrict our considerations to the case a = and rewrite the transport 
parameter Taa in the following form, in which the dependence on temperature and scattering 
phase shift becomes particularly clear [12]: 



Tg{T) = TZ {Caa + (1 - Caa) [sin^ SoVi^ {T) + cos' 6oY(!\T)] } 



(53) 



In (53) Tj^ denotes the normal state limit of Taa- Eq. (53) represents an interpolation procedure 
for the temperature dependence of the transport parameter Taa, which uses the fact, expressed 
in Eq. (51), that at high BQP energies the BQP relaxation time does not need to be evaluated 
sclf-consistently. This gives rise to the definition of a set of generalized Yosida functions of the 
form 




El-^l\DiEp)\^-- y 



(54) 



FS 



which describe the temperature dependence of Taa- On the other hand, an inspection of 
(50) shows, that in the limit Ei^ the elastic relaxation time tends to a finite value in 
the unitary limit, which reflects existence of an impurity band originating from resonant pair- 
breaking processes, described by the t-matrix. This leads to a low-temperature offset in the 
temperature dependence of Taa, which is described by the dimensionless parameter 



Caa — 



2 v"2 



2rN {aj 



'FS 



3/2 



(55) 



FS 



12 




Figure 1 : Temperature dependence of the generalized Yosida functions Yaa {T) , characterizing 
the transport parameters Taa, for the Raman polarizations Big and B23 in the unitary (n = 1) 
and the Born (n = 3) limit, as evaluated from Eq. (54) in the text. The stated power laws fit 
very well at intermediate temperatures. In the low temperature limit, even analytical expressions 
for Y^a\T) are found (c. f. Eqs. (62) and (63)). 
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For a d-wave gap of the form (2) the integrals over the Fermi surface can be performed in Eq. 
(55) in the unitary hmit with the result: 



Cn. 



2 h 



y" 

^0 



TT TnAo Ao 



2 h 

vr TnAo 



K 



E 



(t) 



B 



l9 



(56) 



V" 

^0 

V"2 
^0 



B 



2fl 



B 



l9 



V"2 
1 -I- =2- 



1 _i 4Ao 

2 E" 



B 



2ff 



Here -ftT and £^ refer to the complete elliptic integrals of first and second kind, respectively. 
The description of the transport parameter Taa becomes exact at T = and represents a very 
accurate approximation just below the transition temperature. For intermediate temperatures, 
Eq. (53) represents a physically quite transparent temperature interpolation scheme, which, 

(n) 

however, is meaningful only as long as the generalized Yosida functions Yaa vanish in the low 
temperature limit T ^ 0. We turn now to an evaluation of Eq. (53) in the low temperature 
limit. In the limits of unitary {Sq — 7r/2) and Born (^o 0) scattering, we obtain for not too 
large values of Sq/Aq (i. e. low impurity concentrations): 



limTQ '°=' 



T-»0 



NF{al)FS 



2h 

ttAo 



2h 



In ^ 

^0 



Bis 

B2ff 
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This important result shows an amazing qualitative difference between the Bi^ and the B2g 
polarization: in the Bi^ case, the transport parameter depends on the parameters n,; and 5o, 
characterizing the impurity scattering, whereas in the B2c; case it does not. This behavior of 
the zero temperature transport properties in the case occurs also in the case of electronic 
conductivity [21], the electronic thermal conductivity [22], the (electron-phonon interaction 
induced) sound attenuation [23] and has been termed universal transport [21]. For the case of 
Raman scattering, this result has first been derived for the unitary limit in ref. [24]. 
For the numerical computations it is important to see that in the Born limit (Sq 0) there is 
a simple relation between Sq and the normal state lifetime tn, namely 



Sn = 4Ao exp 



h J 



In the unitary limit, on the other hand, one has to solve the transcendental equation 



II \ 2 



/4Ao 



(58) 



(59) 



^0/ \ ^0 / ^TnAo 
in order to relate Sq to tn. Defining 7 = ^/tnAq, the solution of Eq. (59) can be expressed as 



yll 
^0 



T 

32 / 



(60) 
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with W the Lambert-Ty function, for which an expansion for small arguments reads [35]: 



\W-i{z)\ =ln- +ln ( In 



1 



1 + 



1 

li^ 



+ ... 



(61) 



We consider finally the temperature dependence of Taa in the low temperature limit, where 
an analytical treatment is possible. Using the expansion of D{x) = x[l + i(2/7r)ln{x/4)] for 
X = E'k/Ao <^ 1, one obtains in the Bi^ case 

4 
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Here we have defined the coefficients 
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Figure 3: Normalized imaginary part of the impurity self energy Sg = 9Se(0) in the unitary 
and the Born limit as a function of the normalized scattering rate Tjsf = I/t^v- 



In Fig. 1 we have plotted the generalized Yosida functions Ya" (T), which characterize the 
temperature-dependence of the quasiparticle transport parameters Taa vs. reduced temperature 
T/Tc [c. f. Eq. (52)]. Clearly, the theory predicts a temperature-independent result for Taa 
in the B2g case in the Born scattering limit. This was actually the motivation for Pethick and 
Pines [17] to apply the t-matrix to the description of transport in heavy Fermion superconduc- 
tors, since the experiments showed transport parameters vanishing at low T instead of staying 
constant below Tc. It is quite amazing, that the T-dependence of Yaa\T) is close to the power 
laws T^'^ in the Big- and T^'^ in the B2g predicted by Moreno and Coleman for the ul- 

trasound attenuation [26]. In Fig. 2 we have plotted the normalized parameters CaaTTSS, which are 
proportional to Taa{0), vs. s-wave scattering phase shift cot (5o for Big (op = cos(2(/))) and i?2g 
(op = sin(2</))) symmetry. These parameters characterize the low-T offset in the T-dependence 
of the transport parameters according to Eq. (52). A strong dependence of these offsets on Sq 
and ?iAo/tn is seen in the Bi^ case, whereas in the B2g case Taa{T 0) is independent of 5o 
and /iAq/tn, at least in the limit of low impurity concentrations HAq/t^ — > 0. In Fig. 3 we 
have plotted the dependence of the quantity Sq/Aq, which enters the low-T offset parameter 
Caa in Eq. (53), as a function of the normalized normal state scattering rate Fn/Aq for both 
the unitary and the Born limit of quasiparticle scattering. 
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9 Raman quasiparticle transport and sound attenuation 



In this section wc would like to apply the results obtained in section 8 to the transport param- 
eters connected with the attenuation of ultrasound and Raman scattering, and discuss strong 
similarities in their polarization dependence. The vertex function Op reads in these two cases 
for dimension d: 

{ap = pfvf [(p • q)(p • u) — • u] ; sound attenuation 
7p = mej ■ Mp -65 ; Raman 

Here the unit vectors q and u refer to the sound propagation and polarization directions, re- 
spectively, whereas denote the polarizations of the incoming and reflected photon in an 
electronic Raman scattering process. Without limiting generality, we may assume 

/ cosa \ , . 

q,e/ = . (67) 
\^ sma y 

^ -V / cos (5 

\ smp 

With this 2-d representation, the stress tensor vertex and the Raman vertex assume the form 

ap = ^|cos(a + /?)#?« + sin(a + /3)$p'«} (68) 

7p = cos(a - (3) [const + 7^''^^'] + cos(a + Ph^^'^p^' + sin(a + f^h^^'^p^' 

In (68) the quantities 70 are constants which depend on the band structure and the functions 
$p denote the relevant harmonics of the Fermi surface [27] . In the simplest case of a cylindrical 
Fermi surface one may write 

= cos(4(/)) ; = cos(2(^) ; = sin(20) (69) 

Note that since by definition (a'p)pg = 0, the stress tensor vertex function does not have a 
contribution from the A\g basis function. In other respects, this appears to be the only difference 
between Raman and stress tensor vertex. Note that in particular the dependence of the transport 
parameters on the directions (q, u) and (e/, 65) is the same. With x' = (x + y)/\/2 and 
y' = (x— y) /\/2, this means that the attenuation of transverse sound propagating in the x- (x'-) 
direction corresponds to the xy- (x'y'-) Raman polarization coupling to the 62^- (Big-) basis 
functions. The attenuation of longitudinal sound, on the other hand, propagating in the x- (x'-) 
direction, corresponds to the xx~ (x'x'-) Raman polarization and hence to the combinations 
Aig+Big (Aig+B2g) respectively [27]. We proceed now to write down the transport parameters 
associated with sound attenuation and electronic Raman scattering. In order to be able to 
compare these transport phenomena, we therefore limit the following considerations entirely to 
the Big and B2g basis functions and define generalized Yosida functions (c. f. Eq. (54)) 




A. 



FS 



S = Blg,B2g 

and the dimensionless (offset) parameters (c. f. Eq. (55)) 

c _ " 1 / (•^p)'^" \ . ,_B B mi 
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Let us start with the sound attenuation, which can be described by the wavenumber a{ijj) [36]: 



a(u;) = (72) 
Here ^ is a viscous dissipation parameter, which has the following form 

-.rj.. I {Cs + {1-Cs)yP{T) ■ unitary limit 

Ps{a, P) = cos^(a + P)Ss,Big + sin^(a + P)Ss,B2g 
1 

In (73) r/N denotes the impurity-limited normal state viscosity. The Raman case, on the other 
hand, is characterized by a transport parameter which is reminiscent of the second viscosity, 
discussed in section 6, generalized to include the square of the Raman vertex in the Fermi surface 
average: 

T^^{T) = C^'^{T) (74) 

C77 has the form 

A^)^T^ - ^(3)Nr^^^2^ m/^^ + (^"'^*)^^^^^(^) ' unitary limit 

Hence we have demonstrated the physical similarity of the hydrodynamic transport parameters 
emerging from the electronic Raman effect and the sound attenuation in d=2. 



10 Conclusion 



In summary, we have demonstrated, how the response and transport properties of unconven- 
tional superconductors are linked together in a very simple way, if the long wavelength limit 
q ^ is taken. In this limit the contributions of the normal component (the BQP system) 
and the condensate (the system of Cooper pairs) to the response and transport simply add 
up and allow for a two-fluid description which is valid at arbitrary quasiclassical frequencies. 
The result for the quasiparticle transport parameters Taa^^ and Taa are qualitatively different 
before and after the renormalization of the response with respect to the long-range Coulomb 
interaction. These transport parameters can be represented in a form, which is valid both in 
the limits of weak (Born) and strong (unitary) scattering and which makes their dependence 
on temperature and the impurity scattering parameters {ni,do) particularly clear. In the low T 
limit, the temperature dependence can be evaluated analytically and is found to differ from a 
simple power-law behavior, predicted earlier [29]. In the limit of Born scattering, some compo- 
nents of the quasiparticle transport tensors stay constant in the low T limit and thus exclude 
the possibility of weak scattering as a model for impurity limited transport for example in the 
cuprate superconductors. In the unitary limit, on the other hand, we find offsets in the transport 
parameters, which, in certain cases, do not depend on the impurity scattering parameters. This 
so-called universal transport occurs in both the impurity-limited electronic Raman effect and 
the attenuation of ultrasound. We could finally demonstrate, that for quasi-2-d systems the 
Raman scattering intensity is closely related to the transport of momentum (stress tensor) in 
the BQP system. 
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Appendix: Response and relaxation of densities and currents 



In this appendix we would like to demonstrate that although one has to carefully distinguish the 
response of quasiparticlc densities and currents (their reactive response is indeed different) their 
dissipative response, characterized by their transport parameters, is the same. This justifies the 
line of arguments presented in section 5 of this paper. We now give a more general derivation of 
response functions and transport parameters for a system of Bogoliubov quasiparticles (BQP) 
in (i~wave superconductors. The dynamics of such a BQP system is governed by a scalar kinetic 
equation for the distribution function (5t'k(q,t<j) (c. f. Eq. (21) in the text): 

a;(5z>k — q • Vk/ik = i^Ik (76) 

Here Vk = 9£'k/^9k = (^k/-E'k)vk is the BQP group velocity and the distribution function 

hk = Si^k + VkSEk (77) 

describes the deviation from local equilibrium. In (76) SI^ is the collision integral, which is 
assumed to be limited to the case of purely elastic scattering and to have therefore the following 
(quasiparticle number conserving) form 



Wk = sil+^ + sil'^ 

^k ^kTk /^^ 



(-) _ _K 



(-) 



V 

pa 



Here the index ± denotes the parity of the distribution function with respect to the operation 
k — k. The first part (+) of the collision integral describes the conservation of the BQP 
density, whereas the second part (-) describes the relaxation of the BQP current, here, for 
simplicity, treated within a simple relaxation time approximation. At this stage it is important 
to decompose the BQP energy change SEi^ into even and odd contibutions w. r. t. the operation 
k^ -k: 



5£^k = Ok 



(79) 



In (79) the quantity 6Ca denotes a generalized vector potential. Note that in the case of the 
electromagnetic response (ok = e) one has 6£,a = ^ and = — A/c. The macroscopic BQP 
density is defined through (c. f. Eq. (21) in the text): 

Sn^ = {^P^^'^p) (80) 
Using this definition, one may derive from (76) the following relaxation equation for 6na- 



Here we may identify the BQP current ja in the form 



q-j^ = n«P^wr) (81) 



.(-)\_ / . ip„ u(-) 



j» = {^pi^^ph'p~') = [<^P^JpK~') (82) 
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In the case (ap)FS = 0, the solution of the kinetic equation for the BQP density 5n^ , discussed 
in the text, reads after invoking the effects of the long-range Coulomb interaction (c. f. Eq. 
(24)) 

We wish now to write down the corresponding solution for the generalized current ja . Prom 
(76) one finds for the distribution function /t^"^ 

1 — lUJT^ ^ 1 — lUT^ ^ 

In the long wavelength limit q — this result, inserted into (82) leads to 

= X?».<5C ; Xi(a;) = ^a^ ^vp : vp^^^yp^ (85) 

One immediately recognizes that the response functions for the BQP density Xaa(^) and the 
quasiparticle current Xaa(ci;) differ, besides the tensor structure of the latter, in the usual way in 
their dependence on the dimensionless quantity ujt^. It is instructive to decompose the density 
response function Xaa into its real and imaginary parts, respectively: 

Aaa Aaa ' "Aaa 

^2 - -("pirrl^) w 



X. 



aa ^ \ "p r2 1 I /, „e^2 / " ^-^aa 



In the same way we obtain for xSa- 



aa ' " aa 

2 f, ,^e\2 



^2 - («j|v,:v,J^^ (87) 

It is interesting to formally compare the response of the BQP density and current in the following 
way: The density response can be written in terms of an effective BQP relaxation time Taa in 
the form (c. f. Eq. (33) in the text): 

Sr^Q - _^>3^se ■ - \ ^^^/ (88) 



ation time Tavav 



whereas the current obeys a Drude-like law, involving a corresponding effective current relax- 

-ioJrl^vav / 2 ?P / ^\ Q {^^p^p^^ E^i^^) 

Ja = : -^(«p^p-^yp(Vp-<5Ca)) ; r^vav=J-^ ^.2 ^ \ (89) 

-tu;TS,a„ \ ^p / ((apVp)2^^^\ 



B2 ^ 
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Eqs. (88) and (89) are clearly seen to describe the formal difference between the density and 
current response: whereas the effective relaxation times for the densities and currents differ only 
by the use of different vertices ap ^ flpVp, the response of the current is Drude like, whereas 
the response of the density is not. Nevertheless we have now arrived at a stage where we can 
analyze the connection between the density and current response functions and the transport 
parameters associated with the densities and currents. Defining generalized forces 

fa = -il^S^a ; fa = -i<^Ka (90) 

we may write 

(91) 



= Xg:.(5C + Tg,.f„ 



An inspection of Eqs. (83, 85, 88, 89) shows that while the reactive response of BQP densities and 
currents is qualitatively different, the dissipative response of densities and currents, represented 
by the transport parameters Taa and Taa is similar in structure, if on makes the replacements 
for the vertices Op — OpVp. The transport parameters Taa and Taa read 

Tlio.) = -^ = (ali^^^) (92) 

t2.(^) = 

Hence we have arrived at a unified description of transport phenomena associated with the 
Raman and stress tensor response and have therefore justified the argumentation presented in 
section 5. 
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